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Abstract 

The open electron resonator, described by Duncan et.al 0, is a meso- 
scopic device that has attracted considerable attention due to its remark- 
able behaviour (conductance oscillations), which has been explained by 
detailed theories based on the behaviour of electrons at the top of the 
Fermi sea. In this work, we study the resonator using the simple quan- 
tum quantum electrical circuit approach, developed recently by Li and 
Chen 12]. With this approach, and considering a very simple capacitor- 
like model of the system, we are able to theoretically reproduce the ob- 
served conductance oscillations. A very remarkable feature of the simple 
theory developed here is the fact that the predictions depend mostly on 
very general facts, namely, the discrete nature of electric charge and quan- 
tum mechanics; other detailed features of the systems described enter as 
parameters of the system, such as capacities and inductances. 
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Figure 1: A drawing of the open resonator of Duncan et. al. T, showing the 
reflector (left) and the two point contacts. 

1 Introduction 

The field of mesoscopic physics deals with the frontier between classical and 
quantum physics. However, this frontier is not sharply defined, and every case 
needs special attention. Phenomena like persistent currents [H], Coulomb block- 
ade, magneto-resistance fluctuacions and current magnification ^013 EI] are 
usually studied in this area. The open electron resonator is a device that has 
attracted some interest recently [HI [T] It consists of an open cavity in a 
two-dimensional electron gas, into which electrons can be injected via a quan- 
tum point contact, as described by Duncan et. al and shown in Figure ^ 
Some very interesting properties have been found for this device, such as the 
conductance peaks observed upon changing the cavity length by A//2 (half the 
2D Fermi wavelength). Recently, Duncan et. al 1 have presented a series of 
conductance measurements, showing that the measured conductance through 
the point contact possesses a series of maxima, as a function of the gate voltage 
Vg- These measurements have been performed both at zero magnetic field (B), 
and as a function of it. 

In this work, we present an alternative approach that allows us to explain 
the observed conductance peaks, as a function of the gate voltage Vg, at zero 
magnetic field, B = 0. This approach is related to the analogy between the 
mesoscopic device and a quantum electrical circuit with charge discreteness, 
and it will be discussed in the next section. In sectionOwe study the properties 
of the energy spectrum of our model. 

2 Our model 

We start by observing that the open resonator of reference ^ is similar to 
a three-conductor system, in which one of them, the gate, is kept at voltage 
Vg, while current flows between the other conductors. We complete our model 
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device by assigning to it a self-inductance L. Classically, the circuit may be 
described as a sort charged L — C circuit. We assume then that a current / 
flows between the two charged conductors (with charges Q and — Q); further, 
for simplicity, assume a plane geometry, and compute the classical electrostatic 
energy, obtaining 

U ^C,V^ + ^^+O^QVg, (1) 

where Ci and C2 are capacities of the system, and a is a geometry-related 
number of order unity. Note that this fact will enable us to obtain some rough 
numeric estimates, based on the physical dimensions of the device, which are in 
reasonable agreement with observations, as it will be shown later. 

As it has been proposed by Li and Chen jj], and later by Flores and 
Flores and Utreras |10j , we treat our mesoscopic system as a quantum electrical 
circuit, with quantized charge. The classical Hamiltonian Hd of the (circuit) 
model system, written in terms of the canonically conjugate variables Q and $ 
is 

Hc.-^L^U^- + C,V^ + ^^+aQV,. (2) 

To quantize the system, the variables Q and <& are replaced by the operators 
Q and Further, to recover the quantization of charge within this electrical 
circuit approach, we introduce the replacement |H1 El 

s.f. „(!$), ,3, 

where is the quantum of charge. We remark that, if charge discreteness is 
neglected, the operator <I> may be directly identified with the magnetic flux 
operator, and therefore directly related to the current; however, when one in- 
troduces charge discreteness via the replacement above, the simple relation to 
the current is lost, therefore, after replacement Q the flux operator becomes 
the pseudo-flux. This pseudo-flux operator satisfies the usual commutation 
relation[Q, $] = ifi. Notice that, given the complexity of dealing with an opera- 
tor such as the one above it is simpler to work in the so-called pseudo-flux 
representation, in which the operator $ is replaced by its eigenvalue cf), while the 
charge operator is given by Q = ihd/d(j). In this way, the resulting hamiltonian 
is given by 

H=— sm\q,cj)/2h) + CiKf + ^ + aQVg. (4) 

The hamiltonian operator above constitutes our starting point, and our 
working hypothesis. The parameters of our theory, particularly L and C, are 
related to the geometry of the system, but are hard to compute for a given 
experimental system, therefore, they should be deduced from experimental ob- 
servations. 
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3 Energy spectrum 



To study the properties of our system, define first f3 = (Ci — a^C/2), C = C2, 
and a = aCVg/h, and rewrite our hamiltonian as 

H - + ^ sin'M^h) + ^{^^ + af. (5) 

It is observed that, in this representation, the magnetic energy term acts hke 
a periodic potential in the variable (p, while the electrostatic energy term here 
becomes a kind of 'kinetic' term, i.e., formally, the terms reverse their meaning, 
when one compares with a 'regular' Schrodinger equation. The term /SVg is 
just a constant, and it will be eliminated from the hamiltonian. To make our 
analogies more clear, define a periodic 'potential energy function' V{(l)), with 
period (/iq = 2Trh,/qe = h/q^ 

= ^ sin2(g,0/2a) = \t{1- cos{q,q^/h)} . (6) 
The Schrodinger equation H"^ = E"^ may now be written as 

[VW + ^(*^ + = E^i<l>)- (7) 

The wavefunction ^'(0) satisfies the strong periodicity condition 

+ 00 ) = ^' (8) 

this condition results from the fact that quantization in the charge space implies 
periodicity in the (pseudo) flux space, as shown by Li and Chen 

Observe that the hamiltonian H above commutes with the discrete transla- 
tion operator Tm^g^!{4>) = '^{(p + nupo), for integer m; therefore, for any value 
of a, Bloch's theorem holds, specifically, the energy E{a) should be a periodic 
function of a, labeled by a 'band index' n, further, the parameter a should play 
the role of the wavenumber in Bloch theory; the period for a is Aa = qe/^^■^ and 
therefore, the energies are periodic functions of the gate voltage Vg, with period 
AV^ = q,,/ aC (recall that a is a number that we take as unity from here on). 

To proceed with the numerical calculations, we should get some estimates 
of the parameters of the theory, namely, the values of the capacity C and the 
inductance L. These may be obtained from the knowledge of the linear dimen- 
sions of the system, i.e., of the order of 1/im, and using the classical formula 
C — eolb/d, in which / is taken as the lateral dimension of the cavity, d the dis- 
tance of the contact points and b is the thickness of the 2D electron gas in which 
the device is immersed, now since / ~ d, and 6 is a fraction of a^m, we obtain 
AVg ~ Q.{)2{Volts), a reasonable estimate of the effect. As to the estimate for 
L, it is a parameter related to the size of the circuit that contains the device, 
i.e., a relatively large number L « ^qD. As a measure of this, we have taken 
the ratio of the magnetic potential energy, to the electric energy q^/2C to 

be of order unity, for the numerical calculations described here. 
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Figure 2: Computed energy eigenvalues (first bands) as a function of the 
gate voltage Vg- The ordinate is the energy, on the dimensionless energy scale 
(0o/2i)/(9e/2C)i the abcisa is the dimensionless parameter a(j)o — {2TTC/qe)VG- 
Observe that the ground state 'band' has a very small oscillation amplitude, if 
shown on the same scale than the other curves, it would appear completely flat. 



Notice that we have estimated the capacity of the system using the classical 
formula, but to include small-size corrections it may be more appropriate to 
obtain tha capacity from a microscopic theory, or, at least from Thomas-Fermi 
theory. 

In Appendix we discuss the procedure used to find the energy eigen- 
functions ue{<P) and their corresponding energies for our hamiltonian. The 
solutions are expressed in terms of two normalized and mutually orthogonal 
functions CE{(t>) y Se{4')- The eigenvalue conditions is expressed by ea.f 121(1 . 
which shows explicitly the periodicity of the energies as a function of Vg ■ 

We have computed the eigenfunctions u = UE{(t>) = Wn.Vc by numerically 
solving the eigenvalue equation f Eg. 1211 in Appendix), for the lowest eigenstates, 
as a function of the parameter a = CVg/H. The full wavefunction u(0) is given 
by the linear combination 



UEi<p)^ACEiq^)+BSEi<P), (9) 

with coefficients given by 



B = 



Ac 

V|AcP + |A5|2" 



The results of our energy calculations are displayed in FigureEl which shows 
that the energies are distributed in energy bands with the periodic structure 
already described (period Aa ~ 27r/(/)o). Observe that the first band has very 
small width, which depends on the relation of electric to magnetic energy. 

4 Conductance Oscillations 

From a general point of view, we may expect the same periodicity in the con- 
ductance as that of the energy. If we assume A as the probability per unit of 
time of decaying transitions between two levels {n, n') then the emission power 
is A{E^{n) — E'^{n')), which must be proportional to the conductance G (the 
proportionality constant is the square of the external bias voltage). Since the 
energy spectrum is periodic in AVg — qe/aC, so is the conductance. 

We have shown that the general features of the energy sprectrum do not 
depend on the value of the inductance L, in the sense that, as long as L 7^ 0, 
the energy spectrum will display the Vg periodicity. However, for small values of 
L, this behaviour will not be seen, since the amplitude of the energy oscillations 
will be very small. This is actually observed (see figure for the lowest energy 
'band', even for moderate values of L, in other words, the details of the energy 
spectrum do depend on L, but its general properties, such as the periodicity 
depend only on the capacity of the system, as shown by the experiments. 

We remark that our approach is able to explain the experimentally ob- 
served ^ oscillatory behaviour of G, as a function of Vg, but the linear de- 
pendence of the current on Vg, is not directly explained by our approach, or, to 
put it differently, the present theory does not describe the linear dependence of 
the current on the gate voltage. We observe that, for fixed Vg, we may change 
the capacity of the system, C" = aC2 , by changing one of the linear dimensions 
of the cavity(/, say), while keeping constant both the distance from the contact 
points constant (d), and the other linear dimension (b); it is clearly seen then 
that the current is also periodic in I, i.e., since we may write for C" 

C = eoalb/d, (10) 

in this way, the oscillatory behaviour of the current (and, hence, of the conduc- 
tance too) as a function of the resonator cavity length may also be explained. 

Now we make an estimate of the size of the effect, as suggested by our theory. 
Note that in C above, the area Z x 6 is the product of the linear dimensions of the 
reflector, while d is the distance to the contact points. We see that two of these 
distances are of order 1/im, therefore, C' eq x b (b being the perpendicular 
dimension to the 2D electron gas, of order l/im, then AVg ~ 0.02(T^), a pretty 
good answer, given the rough estimate for the capacity C used here. 
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5 Discussion 



This article offers an alternative point of view on the open electron resonator, 
while recognizing that there exists an accepted explanation of the effect. Nev- 
ertheless, our alternative description is very simple, powerful and it does not 
require elaborate assumptions. It is a new description in terms of elements like 
capacitors and inductances, which is not too widely different from more usual 
treatments for mesoscopic systems |13| . 

It is also worth considering the issue of the spatial dimensionality of the 
theory. A casual glance of our equations gives one the impression that the 
theory is one-dimensional. This is not the case, since the theory is expressed in 
terms of charge-flux variables, and no reference is made to the dimensionality 
of the devices it describes; in fact, such spatial dimensionality only appears 
through the parameters of the theory, such as capacities and inductances, just 
as it happens in classical electrical circuits. It is worth noticing that such theory 
is more correctly described as "topological", and not "geometrical. Therefore, 
a correct use of the theory may help in sheding some light in other mesoscopic 
systems in study today. 
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A Solution of the Eigenvalue equation 

Consider now the (numerical) solution of the resonator equation (jT)), for the 
eigen- function '^{(j)) and eigen-energies, as a function of the parameter Vg (a), 
with the periodic boundary conditions ©. As indicated previously, the wave- 
function ^'(0) may be written as 

This replacement has the property that it eliminates the constant 'a' from 
the Schrodinger equation for the unknown function u{(j)). This wavefunction 
is not periodic, but it satisfies a simple Bloch-type boundary condition (which 
follows from the periodicity of ^f). 

u((/. + (/.o) = e-*'^^°w(</)) (11) 

Now, to actually find numerical solutions of the Schrodinger equation, and 
their corresponding energies, we proceed as follows. For each value of the energy 
E, we find two independent numerical solutions, which we call C{(j)) and S'((/>). 
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For example, for a fundamental period with endpoints at ^ = and (j) = 4>q, 
these functions may be chosen to satisfy 

Ce{<Po/2) = 1 C^(</>o/2)=0 (12) 

5i5(0o/2)=O 5^(</.o/2) = l (13) 

The wavefunction u may now be expressed as the linear combination 

UE{cj>) = A ■ Ce{cI>) + B ■ Se{<P). (14) 
The boundary conditions are now 



ue{<I>o) = e-'"'^°WB(0) 

u'M = e-'»^°<(0), (15) 
which lead to two linear equations for the coefficients A and B, 



AcA + AsB = 

A'cA + A'gB = 0, (16) 

in which the coefficients are 



Ag = e''"l'«/^S{(l)o) -~ e-""^°/^S{0) 

A'c = e""^°/'^C'{<i)o) - e-""''°l'^C'{Q) 

A'g = e''"'"'/^S'{M-e-''"<"'/^S'{0) (17) 

The eigenvalue equation for the energies is given by the vanishing of the 
determinant of the coefficients, 

G{E) = AcA's - AsA'c = (18) 



When, as it happens in our case, the potential is symmetric with respect 
to the midpoint of the interval, the functions C and S have definite parity 
about that point, i.e. C{(j)o) = C(0) and S{(f)f)) = —S{0), the derivatives have 
the opposing parity, C"(^o) = ~C''(0) •S'X^o) = Sl'{Q), therefore, defining 
= a^o/2, we obtain a simplified eigenvalue equation, 

G{E) = sin2(a</)o/2)C(0)S"(0) + cos2(a(/)o/2)S'(0)C"(0) = (19) 

This equation is convenient for numerical work, since the functions C and S 
are easily computed numerically; to find the eigenvalues E, we fix a and compute 
G{E), the zeros correspond to the eigenstatcs. 

For the general case (non symmetric), it is useful to define the coefficients 
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= C(<^o)-C(0) 




= C((/)o) + C(0) 


Si-) 


= S{cf,o)-S{0) 




= S{(f>o) + S{0) 


^(-) 


= C"(</>o)-C"(0) 




= C'{cf>o) + C'{0) 


^(-) 


= S'{ct>o)-S'{0) 




= S'{cj>o) + S'{0). 



With this, the coefficients of the eigenvalue equation become 



The eigenvalue equation is purely real (its imaginary part $J[G'(i<^)] = van- 
ishes identically since the Wronskian -H^[C, S]- of two solutions is a constant): 



G{E) = - ^(-)C('_)) cos2(0)+(5(+)C(V) - sm\e). (21) 
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